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ABSTRACT 



Context. Given the importance of white dwarfs (WDs) in many fields of modem astrophysics, the precise knowledge of the actual 
degree of accuracy of the associated theoretical predictions is a primary task. In the first paper of a series dedicated to the modeling 
of WD structure and evolution we discussed the limits of the available theoretical studies of cooling sequences. 
Aims. In the present work we extend this analysis to isochrones and luminosity functions of WDs belonging to old stellar systems, 
like globular or old disk clusters. The discussion is focused on the most common DA, those with a CO core and an H-rich envelope. 
Methods. We discuss, in particular, the variation of the age derived from the observed WD sequence caused by different assumptions 
about the conductive opacity as well as that induced by changing the carbon abundance in the core. 

Results. The former causes a global uncertainty of the order of 10% and the latter of about 5%. We discuss different choices of the 
initial-to-final mass relation, which induces an uncertainty of 8% on the GC age estimate. 

Key words. Stars: white dwarfs - Stars: luminosity function, mass function 



1. Introduction 

In the last decade, impressive improvement in observational 
techniques made available a significative sample of cold and 
faint white dwarfs (WDs) both in the field (Bergeron, Leggett 
& Ruiz 20011 and in stellar clusters (Paresce, De Marchi & 
Romaniello 1995 Cool, Piotto & King 1996, Renzini et al. ll996[ 
Richer et al. ,1997. Von Hi ppel & Gilmore 2000 Z occali et al. 
|200n H ansen et al. |2002l De Marchi et al. |2004l Monelli et 
al. 120051 ). The growing amount of data prompted a renewed in- 
terest in theoretical studies of WD evolution, particularly for 
the latest stages when the core crystallizes and the stars ap- 
proaches the Debye regime (Benvenuto & Althaus [1999l Hansen 
[19991 Ch abrier et al. l2000l Salaris et al. l2000l Prada Moroni & 
Straniero'2002', hereafter Paper I). 

Concerning the core composition, the theory of stellar evolu- 
tion predicts three kinds of WDs: helium, carbon/oxygen (C-O), 
oxygen/neon/magnesium. The most common are, by far, those 
with a C-O core, which represent the final stage of the evolution 
of single (non-interacting) stars with mass lower than about 7 
Mq (Mup) and a lifetime shorter than the age of the Galaxy. The 
envelope composition is also very important for the evaluation 
of the cooling timescale. Mass loss, particularly during the AGB 
and post- AGB, may substantially erode the envelope. Moreover, 
the H-rich envelope can be substantially reduced in the born- 
again evolution scenario (Iben et al. ,198 3 , Iben & MacDonald 
119951 Herwig I2005I I. Nevertheless, the large majority of the 
observed WDs show an hydrogen-rich atmosphere (DA WDs), 
while only 20% of them are hydrogen-deprived (DB WDs, see 
e.g. Koester & Chanmugam 1983 ). The present study deals with 
the most common WDs, namely the DA type having a C-O core. 
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Here we are interested in the use of WDs as cosmic 
chronometers (Schmidt 1 19591 ) and standard candles (Fusi Pecci 
& Renzini il979i . Due to the intrinsic faintness of WD sequences 
in stellar clusters, only very recently have these tools become 
useful. Their application to the nearest globular clusters (GCs) 
provides age and distance estimates that are independent of those 
obtained by means of the classical methods (those based on the 
main sequence and horizontal branch stars), thus allowing us to 
further constrain the lower limit of the Universe's age (Hansen 
et al. 2002 De Marchi et al. 2004). One should, however, check 
the reliability of the theoretical cooling time predictions before 
adopting WDs as cosmic clocks. Notwithstanding the signifi- 
cant advances in the computation of the WD structure and evo- 
lution, there is still a sizeable uncertainty in the prediction of 
the cooling times of the old WDs (Montgomery et al. 119991 
Salaris et al l20001 Fontaine et al. ,200 li Paper I). At low lumi- 
nosity (logL/LQ 5.5) the differences among the most recent 

models available in the literature reach about 4 Gyr (see e.g. fig. 
1 in Paper I). In Paper I we analyzed some of the main sources 
of uncertainty in the evolution of a 0.6 Mq C-O WD. In partic- 
ular, we discussed uncertainties due to both the WD progenitor 
history and the adopted input physics. 

The evolution of a WD can be roughly described as a cool- 
ing process where the luminosity is supplied by the heat content 
of the internal matter; thermonuclear energy generation is neg- 
ligible, except during the very early phase. Thus, the evolution 
results from the balance between the thermal energy stored in the 
C-O ions constituting the core (about 98% of the WD mass) and 
the energy transport through the He-rich mantel and the H-rich 
envelope, the most opaque region of the star 

In Paper I, we showed that the factors that most severely limit 
the accuracy of theoretical predictions are the assumed conduc- 
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Fig. 1. WD isochrones for ages in the interval 8-14 Gyr in the theoretical plane, i.e. logL/LQ vs. logTe (panel a); V(606) vs. U(336)- 
V(606) (panel b); V(606) vs. B(439)-V(606) (panel c); V(606) vs. V(606)-I(814) (panel d). We adopted the HST/ACS transmission 
curves F336W (broad U), F439W (broad B), F6Q6W (broad V) and F814W (broad /). 



live opacity in the partially degenerate regime, a condition that is 
usually fulfilled in the outer layer of the core and in the He-rich 
mantel, and the amount of carbon and oxygen in the core. At 
logL/LQ ~ -5.5, these two sources produce a total uncertainty 
on the predicted cooling age of about 27%. However, the analy- 
sis reported in Paper I gives only a rough idea of the uncertainty 
in dating GCs with WDs, while a more precise estimate of this 
uncertainty should be based on WD luminosity functions instead 
of tracks of given mass, since this is the best WD chronometer 
For this work we have calculated several grids of models of C-O 
WDs with a hydrogen rich atmosphere, with mass in the range 
0.5-0.9 Mq, under different assumptions of the conductive opac- 
ity and the C-O profile. The related isochrones and luminosity 
functions (LFs) are compared to reveal the main uncertainty fac- 
tors. 



2. Dating WDs in old stellar clusters: the piling up 
in the luminosity function 

For old stellar systems, such as GCs, the best chronometer rely- 
ing on WD evolution is their luminosity function instead of the 
related isochrones. Figure [1] shows WD isochrones for ages in 
the interval 8-14 Gyr in the theoretical plane (logL/LQ vs. logT,,, 
panel a) and in three CM diagrams for different HST/ACS pass- 
bands, F336W (broad U), F439W (broad B), F6Q6W (broad V) 
and F814W (broad /). In the theoretical plane (panel a), there is 
a clear and unambiguous age dependent feature: the blue-turn of 
the WD isochrone, which plays the same role of clock pointer as 
the turn-off point in the main sequence phase. The luminosity of 
the blue-turn depends only on the evolutionary time scales, both 
of the WDs and of their progenitors. On the other hand, on the 
observational CM diagram, the situation is less clear, in fact the 
morphology of the old (>10 Gyr) WD isochrones significantly 
depends on the adopted pass-bands, particularly the luminosity 
of the blue-turn. In the y(606) vs. y(606)-/(814) diagram, com- 
monly used in observations of cool WDs, the blue-turn is almost 
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Fig. 2. Theoretical LFs for 12 (solid line), 13 (dashed line) and 
14 Gyr (long-dashed line). Each magnitude bin is AMv=0.1 
mag. The total number of WDs is 10000. 



insensitive to the age for isochrones older than about 10 Gyr 
This behavior is a consequence of the strong depletion of the IR 
flux in the emerging spectrum of cool (T^ < 5000 K) DA WDs 
due to the collision induced absorption (CIA) of H2-H2 (see e.g. 
Bergeron, Saumon & Wesemael 1995). Since the latter blue shift 
is an atmospheric effect that reflects the departure of the emer- 
gent spectrum from the black body appearance, its extension de- 
pends sensitively on the wavelength of the transmission curves 
of the filter adopted in the observations. The blue shift will be 
much larger in the V-I color than in the B-V, because the CIA is 
the main opacity source in the IR region of the spectrum. Such 
an occurrence explains why in the V vs. V-I plane the blue-hook 
luminosity of old WD isochrones is essentially insensitive to the 
age and the blue-tail of the 14 Gyr WD isochrone is brighter than 
the 12 Gyr one. The degeneracy of the position of the blue-turn 
with age and the necessity to accurately model the CIA opacity 
sources in the WD model atmospheres severely reduce the effec- 
tiveness of the WD isochrones in dating old stellar systems such 
as galactic GCs. A much firmer estimate relies on star counts, 
i.e. luminosity functions, which directly reflect the evolutionary 
time scales. 

Figure|2]shows the WD LFs for 12, 13, 14 Gyr These and the 
following WD LFs have been computed with Monte Carlo sim- 
ulations that distribute 50000 WDs along the related isochrones. 
This figure clearly illustrates the potential of WD LFs in dating 
old stellar systems. Note that the peak of the luminosity function 
shifts by about 0.3 mag per Gyr in the V band, while the main 
sequence turn-off luminosity, the clock classically used to date 
stellar clusters, shifts only by about 0. 1 mag. This implies that 
the WD age estimate is significantly less affected by the uncer- 
tainty on the distance modulus, which represents the main source 
of error in dating GCs with the method based on the turn off lu- 
minosity. Let us briefly recall the main WD evolutionary phases 
in order to explain the principal features in the WD isochrones 
and LFs. The evolution of WDs is initially characterized by a 
rapid drop of the luminosity and of the central temperature. The 
main energy loss characterizing this initial phase is the emission 
of neutrinos. The corresponding zone of the isochrone closely 



Fig. 3. Initial-to-final mass relationship for the four labeled 
metallicities from our stellar models. See text for details. 

follows the track of a WD of fixed mass. Then, as a WD cools, 
the Coulomb interaction between ions becomes progressively 
greater and leads to the crystallization of carbon and eventu- 
ally oxygen (Abrikosov |19601 Kirzhnits 1960; Salpeter[l96lJ- 
Due to the latent heat released by the liquid-solid transition, the 
cooling rate slows temporarily (D'Antona & Mazzitelli 11990! 
Fontaine et al. l2001l l. At the same time, the external convective 
zone penetrates the region where thermal conduction by degen- 
erate electrons is very efficient. Such an occurrence (known as 
convective coupling), initially produces a further slowing down 
of the cooling timescale (see e.g. Chabrier et al. i2000i Fontaine 
et al. 1200 11 1. The steep increase in the WD luminosity function 
is the effect of such a decrease in the cooling rate. Finally, the 
WD enters the Debye regime, where the heat capacity decreases 
as T^, thus depleting the WD of its main energy reservoir As a 
consequence its luminosity quickly drops. 

The morphology of the luminosity function reflects the vari- 
ations of the cooling timescale: the LF is generally sparse at the 
brightest magnitudes and presents a steep rise, followed by a 
sharp cutoff, in the faintest part (see e.g. Fontaine et al. 120011) . 

The best way to date stellar systems through WDs is the anal- 
ysis of their LF. 

3. The model 

The present WD theoretical tracks have been computed with the 
FRANEC (Chiefs & Straniero 1989), a full evolutionaiy code 
that adopts the Henyey method and that is able to model the 
whole evolution of stars from the hydrostatic pre-main sequence 
phase up to the WD stage. 

3.1. Progenitors and starting model 

To provide the starting model for the WD evolution, we started 
from the ZAHB (Zero Age Horizontal Branch) model. The first 
model consists of a helium core (99% of the total mass) and a 
thin hydrogen rich envelope (Y=0.24 Z=10 "^). Then, we fol- 
lowed the evolution at constant mass, through the core He- 
burning. Due to the reduced envelope, at the central helium ex- 
haustion the model rapidly moves towards the blue region of the 
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HR diagram. During this transition the He- and the H-burning 
shells are eventually active. As discussed in Paper I, the model 
settles on the cooling sequence when the masses of both the He- 
rich and the H-rich layers are reduced to certain critical values 
and the thermonuclear burning dies down. 

We are aware that in the real world a DA C-O WD of a 
given mass can be formed through several different evolutionary 
paths, depending on the mass loss rate and on the initial chem- 
ical composition. This was one of the main points discussed in 
Paper I, where we studied three extreme cases, and we showed 
that the internal structure of the resulting WD is poorly affected 
by the different evolutionary histories (see e.g. Fig. 2 and 3 of 
Paper I). Moreover, the computation of the evolution through the 
thermally pulsing asymptotic giant branch (TP-AGB), although 
more realistic, is very time consuming. Thus, the method we 
have adopted represents a useful way to provide realistic WD 
starting models consistent with the evolution of the progenitors. 
In particular the chemical abundance profiles of carbon and oxy- 
gen in the core is largely independent of the evolutionary his- 
tory. Obviously, the computation through the TP-AGB evolution 
should be preferred if the initial, hot phase of WD evolution is 
analyzed. In such an early phase of the cooling the outer layers 
still bear the record of the complex chemical abundance pro- 
file left by TP-AGB evolution, which is subsequently modified 
by diffusion. The metallicity of the progenitors determines the 
maximum value of the envelope mass, but this quantity does not 
significantly affect the cooling timescale (see paper I). The only 
important quantity is the relation between the progenitor mass 
and the WD mass (initial-to-final mass relation) which is needed 
to correctly evaluate the time spent before the cooling sequence 
and, in turn, correctly locate a WD with a certain mass on the 
HR diagram (see 2.3). We stop the progenitor evolution when 
the newly formed WD reaches log L/L© ~ 0. At this point, the 
helium and hydrogen in the envelope are fully sorted, in order 
to reproduce the observed chemical stratification of the atmo- 
spheres of DA WDs. Several studies have indeed shown that in 
such a compact stars, where the surface gravity g is of the order 
of ~ 10^ - 10^ cm sec"^, the gravitational settling is efficient 
enough to account for the monoelemental feature of WD spec- 
tra (F ontaine & Michaud 1979 Muchmore 1984 Paquette et al. 
11986) . The model obtained with this method is our starting point 
for the computation of the cooling sequence. 

3.2. Cooling tracks 

Let us recall the major features of the present models; for a de- 
tailed description see Paper I. We put considerable effort into 
updating the main input physics, namely the equation of state 
(EOS) for a high-density plasma, the radiative and conductive 
opacity and the model atmospheres. We updated and extended 
the EOS for a fully ionized plasma originally described by 
Straniero ( I1988I I. In particular, we have adopted the most re- 
cent results for the electrostatic corrections both in the liquid 
and in the solid phase (see Paper I for references). Concerning 
the Rosseland radiative opacity, we adopted the OPAL (Iglesias 
& Rogers 1996) for high temperature (log T[K] > 4.0) and the 
Alexander & Ferguson (1994) molecular opacities for the low 
temperatures (log T [K] < 4.0). 

As we have shown in Paper I, for the cooling evolution of 
old WDs, the treatment of electron conduction in the partially 
degenerate regimes typical of the outer layers is complex. We 
have computed sets of theoretical cooling sequences and the re- 
lated isochrones adopting the conductive opacity by Hubbard & 
Lampe (119691 hereafter HL69), by Itoh and coworkers (Itoh et 



al. [19831 Mitake et al.IMland Itoh et al. [19931 hereaft er 193) , 
Potekhin ([T9991 hereafter P99, see also Potekhin et al. [T999l l. 
and combinations. The description of the different assumptions 
and the related effects on the cooling sequences of WDs is given 
in section 5. 

The computations of reliable cooling models rely on the 
choice of appropriate boundary conditions, namely the pres- 
sure at the base of the photosphere. As shown by Fontaine et 
al. (11974) and then confirmed by successive studies (see e.g. 
Tassoul, Fontaine & Winget 1990' Hansen 1 1999) it is important 
to use a detailed model atmosphere to fix the surface boundary 
conditions. It is well known that one of the most important events 
in the evolution of old WDs, the coupling between the outer con- 
vective envelope and the isothermal degenerate core (see e.g. 
Chabrier et al. 2000, Fontaine et al. 2001) is quite sensitive to the 
details of atmosphere stratification, in particular a gray approx- 
imation overestimate the inward penetration of the convection. 
The quoted coupling initially produces a sizeable slow down in 
the cooling evolution followed by a rapid drop in the luminosity. 
Since the occurrence of the convective coupling, the core, which 
is essentially the energy reservoir of the star, and the atmosphere, 
where the energy is radiated in space, are directly coupled. The 
pressure at the outer boundary, namely the layer where the op- 
tical depth T - 50, has been obtained by interpolating on a grid 
of model atmospheres (Bergeron 2001 ) which cover a range of 
surface gravities between lO^'^-lO^ cm/s^ and effective temper- 
atures between 1500-100000 K. A detailed treatment of the at- 
mospheric structure is fundamental also for the prediction of the 
observable properties, the emergent spectrum and the photomet- 
ric colors. As shown by Bergeron, Saumon & Wesemael ( 19951 1, 
when a WD with a pure hydrogen atmosphere cools down to 
T(. ~ 5000 K and the molecular recombination begins, the col- 
lisionally induced absorption (CIA) of H2-H2 progressively be- 
comes the most important opacity source for the out-flowing in- 
frared photons, so that the emergent spectrum departs from the 
blackbody appearance and the colors are blue shifted (see e.g. 
also Hansen [T999l ). 

The present models take into account the liquid-solid phase 
transition and the related energy release (see e.g. Paper I) but ne- 
glect the effect of the phase separation during the crystallization 
of the core, i.e. the chemical redistribution of carbon and oxygen 
due to the shape of the phase diagram for this binary mixture 
(Stevenson |1980l l. We are interested in an evaluation of the rela- 
tive effects on WD evolution of different physical aspects rather 
than to provide a reference standard model. Since the effect on 
WD evolution of chemical fragmentation has been already dis- 
cussed several times in literature (see e.g. Segretain et al. 1994) 
Salaris et al. 11997 . Montgomery et al. |1999t Isern et al. 1 19971 
2000) we do not further investigate this issue. 

3.3. Initial-to-final mass relation and initial mass function 

In order to compute isochrones, the age of the progenitor at the 
entrance of the cooling sequence must be added to the cooling 
time. This quantity depends on the mass and the chemical com- 
position of the progenitor The relationship between the main 
sequence mass and the final WD mass is still uncertain; it is sen- 
sitive to many poorly known phenomena, like the extention of 
the H-burning convective core, which affects the stellar lifetime, 
and the efficiency of the mass loss mechanism, particularly dur- 
ing the TP-AGB phase when the star loses the greatest fraction 
of its mass, which affects the final mass. 

Semiempirical initial-to-final mass relations are only avail- 
able for the chemical composition of the solar neighborhood and 
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Table 1. Initial-to-final mass relationship for the four labeled metallicities from our stellar models. See text for details. Masses are 
in Mq 



Mms 0.8 0.9 1.0 1.1 1.2 1.5 1.8 2.0 2.3 2.5 2.7 3.0 4.0 5.0 6.0 7.0 8.0 

Z=0.02 

M„D 0.55 0.56 0.565 0.57 0.58 0.59 0.60 0.61 0.62 0.63 0.64 0.67 0.81 0.91 0.985 1.06 1.14 

Z=0.006 

MwD 0.55 0.56 0.565 0.57 0.58 0.59 0.60 0.61 0.62 0.64 0.67 0.80 0.895 0.96 1.04 1.14 1.21 

Z=0.001 

M„D 0.56 0.57 0.575 0.58 0.59 0.61 0.62 0.64 0.68 0.72 0.76 0.82 0.90 0.98 1.05 1.15 1.22 

Z=0.0001 

MwD 0.575 0.59 0.595 0.605 0.615 0.625 0.63 0.66 0.72 0.77 0.815 0.86 0.92 0.99 1.09 1.18 1.26 



for the Magellanic Clouds (see e.g. Weidemann & Koester [l983l 
Weidemann [T9871 |2000l Herwig [T997t . For the typical compo- 
sition of the halo, the long distance and the large age limits a 
direct determination of the initial-to-final mass relation. Thus, in 
spite of the great uncertainty, we are forced to use theoretical re- 
lations. Figure[3]shows the initial-to-final mass relationship from 
the homogeneous set of stellar models of low and intermediate 
mass stars computed by Straniero et al. (1997) and Dominguez 
et al. (|T999), with the same stellar evolutionary code used to 
provide the present set of WD models. The mass loss during the 
AGB phase has been obtained by best fitting the period-mass 
loss diagrams obtained by Whitelock et al. (I2003I I. In stars with 
initial mass lower than 3 Mq, an additional pre-AGB mass loss 
has been considered based on the Reimers formula with the pa- 
rameter rj = 0.4. Table [T] reports the initial-to-final mass rela- 
tions shown in figure |3] Three distinct regions with different 
slopes may be distinguished. Concerning the solar metallicity, 
these three regions are: M,- < 3 Mq, 3 Mq < M, < 4.5 Mq and M, 
> 4.5 Mq, respectively. At lower metallicity the limit in mass of 
these three regions decreases, but the general features are con- 
served. For all the three regions the adopted mass loss prescrip- 
tion during the TP- AGB plays a crucial role in the determination 
of the final WD mass. Stars belonging to the first interval have 
a radiative core during the main sequence, or develop a small 
central convective zone (<< 0.5 M©). Thus the mass of the hy- 
drogen exhausted core (Mh) at the first thermal pulse depends on 
the efficiency of the H-burning shell operating during the RGB 
phase and the following He-burning phase, and on the corre- 
sponding evolutionary time scales. These parameters are very 
similar for all these low mass stars so that their Mh at the first 
thermal pulse is practically the same (~ 0.55 Mq, see e.g. Fig. 
15 of Dominguez et al. |1999l l. Then, the most important phe- 
nomenon that influences the final WD mass of these stars is the 
mass loss, taking place during the RGB and AGB phases, which 
erodes the envelope and limits the further growth of Mh during 
the thermally pulsing AGB. At variance, stars belonging to the 
second zone develop a sizeable convective core when they burn 
the central H. Thus, Mh at the first thermal pulse is mainly de- 



termined by the maximum extension attained by the convective 
core during the main sequence. In this case, the pre-AGB mass 
loss have a negligible effect on the final mass, while that in the 
TP- AGB phase is, once again, crucial. Similar considerations are 
valid for the stars belonging to the third region, but, in this case, 
the growth of the core is also limited by the occurrence of the 
II dredge-up, during the early- AGB phase. During the TP- AGB 
phase the recursive occurrence of a deep III dredge-up brings 
inward fresh hydrogen and puts backward the outer edge of the 
hydrogen exhausted core, thus preventing a substantial increase 
of Mh- For this reason, the final mass is only slightly larger than 
Mh at the beginning of the TP-AGB phase. As shown in Fig. 20 
of Dominguez et al. (1999 ) there is a general agreement among 
the available predictions for Mh at the first thermal pulse as a 
function of the initial mass. The most important consequence of 
the assumption of a sizeable convective core overshoot is the de- 
crease of the maximum mass able to produce a C-O WD (M,,^, 
Bertelli, Bressan & Chiosi 1985 ). In order to compute WD LFs 
and synthetic color-magnitude diagrams (CMDs) we also need 
the initial mass function (IMF) for the progenitor stars. In the fol- 
lowing, if not explicitely stated, we have adopted a single power 
law IMF (dN/dM, oc Mr"), with a Salpeter exponent (a=2.3). 
The effect of varying both the slope of the IMF and the initial- 
to-final mass relationship will be discussed in section 7. 



4. Carbon and oxygen abundances in the core 

The thermal energy of a cool WD is almost entirely stored in 
the carbon and oxygen ions of the core. The cooling rate de- 
pends sensitively on the relative amount of carbon and oxygen. 
Carbon has an higher specific heat than oxygen, simply because 
the number of ions in a gram of C is larger. Thus, the cooling 
rate depends sensitively on the C/O ratio: the higher the carbon 
abundance in the core, the slower the cooling of the WD (Salaris 
et al. 1997; Chabrier et al. 120001 Paper I). In addition, since car- 
bon crystallizes later than oxygen, the release of the latent heat 
and the onset of convective coupling are delayed when carbon is 
overabundant. 
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Fig. 4. Theoretical LFs for 12 (thin lines) and 14 Gyr (thick lines) 
for WDs with low (solid lines) and high (dashed lines) C abun- 
dance. 



The chemical abundance profiles in the core of a C-O WD 
are the result of the helium burning that occurred in the core 
and in the shell of the progenitor star As is well known, initially 
the 3a reaction produces carbon and afterwards the '^C(q', yY^O 
synthesizes oxygen. The final chemical stratification in the core 
of a WD depends on the relative efficiency of these two reac- 
tions. In particular, the reaction rate of the '^C(a, y)'^0 at an 
energy relevant for the stellar He-burning (about 300 KeV) is 
only known within a factor of 2 (see e.g. Buchman'1996' Kunz 
et al. 120021 ). The treatment of the convective mixing during the 
core He-burning phase also affects the predicted C/O ratio (see 
Imbriani et al. 2001 , Straniero et aL lIOOgT l. 

In Paper I we analyzed the changes of the cooling timescale 
of a 0.6 Mo CO WD caused by the variation of the ^^C(a, yy^O 
reaction rate, within the present error bars, and by different pre- 
scriptions for the core convection. 

In this section we extend that analysis to masses from 0.5 
to 0.9 Mq and calculate the corresponding isochrones and LFs. 
Then, two sets of WD cooling tracks were obtained for two dif- 
ferent C-O chemical stratifications in the core, the low carbon 
abundance model with a central carbon mass fraction of about 
0.2 and the high ones with 0.50, a range that roughly represents 
the present uncertainty in the theoretical prediction. The former 
value appears in better agreement with the abundance derived 
from the power spectra of variable WDs (Metcalfe, Winget & 
Charbonneau l2001l l. although the feasibility of this kind of mea- 
surement has been questioned (Fontaine & Brassard , 2002i) . 

Figure |4] shows the 12 (thin lines) and 14 Gyr (thick lines) 
LFs for WDs with low (solid lines) and high (dashed lines) C 
abundance. As expected, due to the lower heat capacity of the 
carbon-poor material, the cooling time is shorter in these models 
than the cooling time of the carbon-rich models. As a conse- 
quence, the WD LF peak occurs at lower luminosity. The effect 
of the present uncertainty on the chemical profile in the core is 
sizeable. 

Taking into account that the position of the WD LF peak is 
a very sensitive function of the age (AM,yAAge ~0.3 mag/Gyr 
for GCs ages), the poor determination of the C abundance in the 
core translates directly into an uncertainty of about 0.6 Gyr in 
the inferred GC age. 
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Fig. 5. Theoretical LFs for 12 (thin lines) and 14 Gyr (thick 
lines) for WD models computed adopting the conductive opac- 
ity of Itoh and coworkers (193) in the fully degenerate regime 
(e = T/r/r < 0.1) and Hubbard & Lampe (HL69) in the partially 
degenerate one (0 > 0.1) (I93-I-HL69, sohd hnes) and 193 in the 
whole structure (193, dashed lines). 
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Fig. 6. Theoretical LFs for 12 (thin lines) and 14 Gyr (thick lines) 
for WD models computed adopting the conductive opacity of 
Potekhin (P99) and HL69 in the regions where, respectively, < 
0.1 and > 1, and a linear interpolation in the transition zone 
(P99-I-HL69 smooth, solid lines) and the ones computed adopt- 
ing the P99 conductive opacity in the whole structure (dashed 
lines). 



5. Conductive opacity 

The electrons in the core of a WD are fully degenerate. In these 
conditions the energy transfer is largely dominated by electron 
conduction whose efficiency is high enough to keep the core al- 
most isothermal. On the contrary, in the mantel (He-rich) and in 
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the envelope (H-rich) the electrons are only partially or not de- 
generate, so that thermal conduction is less efficient. In bright 
WDs, the thin envelope is the most opaque region of the WD, a 
kind of insulating layer which regulates the temperature decrease 
of the core. The progressive development of the external con- 
vection penetrating the partially degenerate region in cool WDs 
causes the already mentioned convective coupling. Therefore, it 
is of primary importance for the description of the cooling pro- 
cess to understand the detailed treatment of the energy trans- 
port through the external layers (see e.g. D' Antona & Mazzitelli 
[19901 Fontaine et al. l200TT l. 

In our knowledge, there are only two published conductive 
opacities suitable for partially degenerate regimes, the pioneer- 
ing study by Hubbard & Lampe (119691 hereafter HL69) and 
Potekhin ( [19991 hereafter P99, see also Potekhin et al. [19991 ) [^. 
Regarding the fully degenerate regime, in addition to the just 
mentioned works, there are also the fundamental contributions 
by Itoh a nd cow orkers (Itoh et al. il983> Mi take et al. i 1984, and 
Itoh et al. [19931 hereafter 193). 

As discussed in detail in Paper I, the HL69, although nom- 
inally valid, should not be adopted in regions with a high 
Coulomb coupling parameter (F > 10, where F=(Ze)^/kBTa and 
a is the interionic distance) due to an outdated treatment of the 
liquid-solid transition. On the other hand, the 193 computations 
are strictly valid only in the completely degenerate regimes (0= 
Tf/T < 0.1, where Tf is the electronic Fermi temperature), in 
fact they underestimate the contribution of the electron-electron 
interactions which is not negligible in the partially degenerate 
regimes. The only conductive opacity computations suitable for 
the whole WD structure and evolution are those by P99. 

In paper I we showed the dramatic effect on the evolution of 
a 0.6 M0 WD of different assumptions for the conductive opacity 
in the thin envelope. At log L/Lq= -5.5 the adoption of the 193 
in the whole WD structure leads to an age 17% younger than 
that obtained by using HL69 in the partially degenerate regime. 
In this section we extended that analysis to evaluate the effect 
on the WD isochrones and LFs of different assumptions on the 
conductive opacity. 

We have computed five sets of cooling tracks with masses 
in the range 0.5-0.9 and the related set of isochrones and 
LFs. Figure |5] shows the theoretical LFs for 12 (thin lines) and 
14 Gyr (thick lines) for WD models computed adopting the con- 
ductive opacity of Itoh and coworkers (193) in the fully degen- 
erate regime (6 < 0.1) and Hubbard & Lampe (HL69) in the 
partially degenerate one (6 > 0.1) (I93H-HL69, solid lines) and 
193 in the whole structure (193, dashed lines). Notwithstanding 
that the two pairs of sets differ only in the conductive efficiency 
adopted in the partially degenerate regime, roughly correspond- 
ing to the thin external layer, whose mass is less than 1 % of the 
WD mass, the effect is very large causing a shift in the position 
of the LF peak of about 0.3 mag for 12 Gyr and 0.5 mag for 14 
Gyr, which means roughly a difference in age of 1 Gyr and 1.6 
Gyr, respectively. 

Itoh suggested to use a different match between the 193 and 
the HL69 conductive opacity (private communication). Thus, we 
have calculated an additional set of cooling sequences, here re- 
ferred to as I93-HHL69 smooth, where 193 is used for 9 < 0.1, 
HL69 for > I and a linear interpolation is used between the 
two for 0. 1 < < 1 . The variation of the WD LF peak is less 
than 0.1 mag. In addition, we have computed cooling models 
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' We have derived conductive opacity tables by means of the 
Fortran codes provided by A.Y. Potekhin and available at the web site 
( ,http://www.ioffe.rssi.ru/astro/conduct/conduct.htmlj . 



Fig. 7. Theoretical LFs of 14 Gyr for the four labeled metallici- 
ties. 



adopting the conductive opacity of Potekhin (P99) and HL69 in 
the regions where, respectively, 6 < Q.l and 6 > 1, and a linear 
interpolation in the transition zone (here referred as P99-I-HL69 
smooth). To compute the conductive opacity of Potekhin we use 
Fortran codes provided by A.Y. Potekhin. The differences are 
small, as expected from a previous analysis presented in Paper I, 
because the P99 predictions are in good agreement with the 193 
ones where this latter is nominally valid. 

Figure|6] shows the LFs for 12 (thin lines) and 14 Gyr (thick 
lines) for the just presented "P99-I-HL69 smooth" WD mod- 
els (solid lines) and the ones computed adopting the P99 con- 
ductive opacity in the whole structure (dashed lines). The dis- 
crepancy between the two sets of models is quite large, for 
the presented ages the shift of the WD LF peak is about 0.4 
mag, roughly corresponding to 1.3 Gyr As in the case shown 
in Fig. |5] this result proves again the extreme sensitivity of 
the computed cooling times and thus of the theoretical WD 
LFs to the treatment of conductive opacity in the very thin 
(Mf„i,/Miyz) <0.01) envelope, where electrons are only partially 
degenerate. To our knowledge, Potekhin provides the only ho- 
mogeneous conductive opacity computations suitable for the de- 
scription of the whole WD structure available in the literature. 
On the other hand, these results for nondegenerate plasmas, as 
stated by Potekhin himself, "are based on a continuation from 
the degenerate domain (using Fermi-Dirac averaging) and can be 
considered as order-of-magnitude estimates." (see e.g. descrip- 
tion at http://www.ioffe.rssi.ru/astro/conduct/conduct.html). To 
provide firmer predictions of WD cooling times we need more 
precise treatments of electron conduction for partially degener- 
ate regimes. 

At present the choice of the prescription adopted for the con- 
ductive opacity affects the age estimate of GCs at the level of 
10%. 



6. Metal abundance 

In Paper I we computed the evolution of two 0.6 M0 WD models 
produced by progenitor stars with very different metallicity, re- 
spectively Z- 0.0001 and 0.02. The difference between the two 
cooling times are almost negligible (see e.g. Fig. 5 in Paper I), 
about 2% at the faint end of our computations (log L/Lq =-5.5). 
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Fig. 8. Initial-to-final mass relations for the solar chemical com- 
position from: Weidemann ( 120001 filled circles); Girardi et al. 
( 12000 1 dashed line); our relation (solid Une); Hurley et al. (120001 
short dashed line). 



Fig. 9. Theoretical WD LFs for 12 Gyr and different initial- to- 
final mass relations for the solar chemical composition from: 
Weidemann (2000 filled circles); Girardi et al. (2000 dashed 
line); our relation (solid line); Hurley et al. (120001 short dashed 
line). 



The larger the metallicity, the thinner the external H-layer, 
because the H burning efficiency increases with a larger amount 
of CNO. Due to the compact nature of these stellar remnants 
the gravitational settling is so efficient as to produce an envelope 
constituted by two layers of almost pure elements, i.e. a buffer 
of He surrounded by a very thin H surface envelope. 

On the contrary, the WD isochrones and LFs are consider- 
ably affected by metallicity variations because, as is well known, 
the pre-WD evolutionary time scale is a sensitive function of 
the original chemical abundance. In addition, the initial-to-final 
mass relation changes with metallicity. 

In order to estimate the effect of the assumed chemical abun- 
dance, we computed four sets of isochrones and the related 
LFs for metallicities covering a range suitable for population I 
and II stars, i.e. Z=0.0001, 0.001, 0.006 and 0.02. As already 
mentioned, we adopted the pre-WD evolutionary times and the 
initial-to-final mass relationship (see e.g. Fig. [3]) from the homo- 
geneous set of stellar models of low and intermediate mass stars 
by Dominguez et al. ( 1999 ) and Straniero et al. ( 1997). 

Figure|2]shows the theoretical WD LFs for an age of 14 Gyr 
and the four labeled metallicities. The position of the peak of 
the WD LF is a sensitive function of the metal content of the 
progenitor stars. Between Z=0.0001 and Z=0.001, there is a dif- 
ference in the visual magnitude of the peak of about 0.4 mag, 
the same between Z=0.001 and Z=0.006. Thus, in order to use 
WDs to date GCs, one should compute theoretical isochrones 
and LFs with the cluster metallicity. For galactic GCs the typi- 
cal error on [Fe/H] is of the order of 0.1 dex. This means that the 
typical uncertainty on the metallicity measurements only slightly 
affects the estimates of the GC ages based on the WD cooling se- 
quences, provided that WD isochrones and LFs have been com- 
puted with the nominal value of the cluster metal abundance. 

7. Initial-to-final mass relation and IMF 

Figure [8] shows some of the most recent initial-to-final mass re- 
lations for the chemical composition typical of the solar neigh- 
borhood. Weidemann (|2000 . filled circles) derived the WD mass 



by comparing theoretical WD evolutionary models with the sur- 
face gravities and effective temperatures inferred from spec- 
troscopy. Then the initial mass is obtained by comparing the 
cluster turn off luminosity with stellar models of different masses 
(see e.g. Weidemann & Koester ll983l ); thus also the semiempir- 
ical method depends on some theoretical assumptions, like the 
convective core overshoot. As mentioned above and as can be 
seen in Fig. [8] the initial-to-final mass relation is still quite un- 
certain, even in the case of typical solar chemical composition. 

Girardi et al. (120001 ) rely on synthetic AGB models. This 
means that the evolution of the AGB stellar properties (the core 
mass, the total mass, the luminosity or the effective temperature) 
is obtained by semiempirical prescriptions, instead of solving 
the stellar structure equations (see e.g. Groenewegen & de Jong 
IT9931 Bertelh et al. 11^941 Marig o et al . [19961 IT^ Girardi & 
Bertelli '1998'). The Hurley et al. (2000) relation has been com- 
puted with the SSE package, kindly provided by Hurley, which 
uses analytic formulae derived by fitting the stellar models com- 
puted by Pols et al. (il995i Our prescription for the initial-to- 
final mass relation is based on full stellar evolution computa- 
tions, from the pre-MS up to the AGB tip (Straniero et al. 1 19971 
Dominguez et al. 1999 Straniero et al. 2005). The pre-AGB 
mass loss is based on the Reimers' formula (77=0.4), and an em- 
pirical calibration of the mass loss-period relation is used for the 
AGB evolution (see Straniero et al. 12005 1 for details). Thus, the 
three curves reported in figure [8] are representative of the differ- 
ent methods commonly used to derive the initial-to-final mass 
relation. 

Figure |9| shows the WD LFs related to the initial-to-final 
mass relations described in Fig. [8] In all cases we have adopted 
the same set of cooling evolutionary models and the same slope 
of the IMF. We derived the pre-WD ages from the homoge- 
neous set of stellar models of low and intermediate mass stars by 
Straniero et al. (1997) and Dominguez et al. ( 1999), except for 
the Girardi et al. (2000) and the Hurley et al. (2000) relations, for 
which we used the life times provided by their own stellar mod- 
els. The substitution of our relation by the Weidemann ( 12000! ) 
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Fig. 10. Theoretical WD LFs for 12 Gyr and solar metallicity. 
Three different slopes of the IMF (oc M"") have been adopted: 
a-2.3 (solid line); a- 3 (dashed line); a- 1 (dot-dashed line). 

one produces a shift of the peak of about 0. 1 mag, roughly cor- 
responding to 0.3 Gyr The same shift, but in the opposite di- 
rection, is produced by substituting our relation with the Hurley 
et al. ( 120001) one. A larger effect is determined by changing our 
relation with the Girardi et al. (2000) one; in this case the peak 
shifts by about 0.3 mag, which translates into an age error of 
1 Gyr. To illustrate the effect of a different mass distribution of 
the progenitor stars, we have computed WD LFs adopting a sin- 
gle power law IMF (dN/dM, ocMj") varying the slope within a 
wide range, i.e. a= 1, 2.3 and 3. Figure [10] shows the results; as 
expected the variation of the slope of the IMF affects only the 
height of the peak, but leaves its position practically unchanged. 

8. Conclusions 

WD cosmochronology is becoming a feasible tool for dating 
stellar clusters. In this paper we have analyzed the accuracy of 
the available calibration of the WDs age, by computing various 
sets of WD isochrones and LFs under different prescriptions for 
the theoretical ingredients. We discussed the effect on the WD 
LFs, and thus on the globular cluster age derived with these com- 
pact objects, of the assumed conductive opacity, C-O abundance 
in the core, and initial-to-final mass relation and IMF. 

The main effect is caused by the adopted conductive opacity, 
in particular in the partially degenerate regime characteristic of 
the thin outer layer. For the range of ages typical of galactic GCs 
the conductive opacity affects the dating by the order of 10%. 

Concerning the C-O amount in the core, the uncertainty in 
the age is less than 5%. Moreover, this contribution will further 
decrease in the near future thanks to experiments to measure the 
^^C(a, y)'*(9 reaction rate at the energy of interest. 

The IMF does not affect the position of the peak and conse- 
quently the estimated ages. At variance, the initial-to-final mass 
relation plays an important role in dating globular clusters be- 
cause it significantly affects the shape of the WD LF, in particu- 
lar the position of the peak. The uncertainty on the age can reach 
8%. 

Finally, given the sensitive dependence of the evolution of 
the WD progenitors on the metal content, in order to use WDs 



to date GCs, one has to compute WD isochrones and LFs for the 
metallicity of the cluster. In fact, between Z=0.0001 and 0.001 
and between Z=0.001 and 0.006, the difference in the inferred 
age is, for the range of interest, of the order of 10%. 
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